This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 


BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 

• BLACK BORDERS 

• TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SKEWED/SLANTED IMAGES 

• COLORED PHOTOS 

• BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 


IMAGES ARE BEST AVAILABLE COPY. 


As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 


(19) 



(12) 


apdiscfies Patentamt 
European Patent Office 
Office europeen des brevets (1 1 ) EP 0 751 482 A1 

DEMANDE DE BREVET EUROPEEN 


(43) Date de publication: 

02.01.1997 Bulletin 1997/01 

(21 ) Numdro de depdt: 96201 774.5 

(22) Datede depdt: 26.06.1996 


(51) Inta 6 : G06T5/00 


(84) Etats contractants designes: 

. PHILIPS ELECTRONICS N.V. 

DE FR GB NL 

5621 BA Eindhoven (NL) 


Etats contractants designed: 

(30) Priority: 30.06.1995 FR 9507888 

DE GB NL 

(71) Demandeurs: 

(72) Inventeur: Makram-Ebeid, Sherrf 

• LABORATOIRES D'ELECTRONIQU E PHIUPS 

75008 Paris (FR) 

S.A.S. 


94450 Limeil-Brevannes (FR) 

(74) Mandataire: Lottin, Claudine 

Etats contractants designes: 

Soctete Civile S.P.I.D. 

FR 

156, Boulevard Haussmann 


75008 Paris (FR) 


(54) Procede de f iltrage temporel du bruit dans une image d'une sequence d'images numeriques 
et dispositif mettant en oeuvre ce procede) 


(57) Un procede de f iltrage d'une image d'une 
sequence d'images numerisees bruitees comprend 
revaluation d'une equation ^integration fournissant 
I'intensite filtr^e (P t c ) d'un pixel donne [A t (x,y)] dans 
llmage du present (J^ en evaluant ta somme de 
I'intensite bruitee (l t p ) du present (t) et du produit d'un 
facteur de gain (K t c ) par la difference de ladite intensity 
bruitee du present (l t p ) et de I'intensite filtree causale 
(P £ ). Dans cette evaluation, le facteur de gain (K t c ) est 
une fonction du facteur de gain causal (K^) et d'un 
coefficient de continuity (a *p) mesurant la probability de 


continuity d'intensite entre I'intensite bruitee du present 
(1^ et I'intensite filtree causale (P m ). 

Un disposrtif pour mettre en oeuvre ce procede 
comprend un sous-ensemble recevant & son entree 
(100) Techantillon bruite du present et incluant une 
m§moire (10) fournissant une intensity filtree causale 
(P t ?i ), une memoire (17) fournissant un facteur de gain 
causal (K t ?i ) et des moyens de calculs de ['equation 
d'integration. 

Application : Imagerie medicale en fluoroscopie 
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Description 

FIELD OF THE INVENTION 

5 (/invention concerr;e un precede de filtrage temporel du bruit dans upe image d'une sequence dlmages numeri- 

ques. ainsi qu'un dispositif pour mettre en oeuvre ce precede. 

Llnventiontrouve particulierement son application dans le trait ement dlmages medicales formees dans un mode 
f luoroscopique de rayons X, au moyen d'un systerne ou le faisceau de rayons X a une faibie intensite, fournissant sous 
forme de sequence des images bruitees et peu contrastees qui doiverrt §tre filtrees afin d'en extraire le bnit sans en 

to deteriorer les details. 

Le filtrage temporel consiste a lisser un signal monodimensionnel appeie signal tempore!, constitue par les vaieurs 
dlntensite d'un pixel ayant une localisation donnee dans les images de la sequence, en fonction du temps. 

BACKGROUND OF THE INVENTION 

75 

Un procede de filtrage temporel est deja connu par la publication de R.E.KALMAN, intitulee "A new approach to 
linear filtering and prediction problems" dans "Transactions of the ASME, Journal of Basic Engineering, Series 82D, pp. 
35-45, 1960. 

Le filtrage de Kalman est defini par une equation recursive donnant I'intensite filtrees d'un pixel courant d'une 

20 image de fa sequence en fonction d'hypotheses faites A PRIORI, en fonction de I'intensite du pixel de meme localisation 
dans Ilmage precedente de la sequence et en fonction d'un facteur appeie gain de Kalman. 

Cette equation peut conduire a deux algorithmes recursifs. Un probleme est que, des quli se produit un petit mou- 
vement entre deux images, ce mouvement donne lieu a un front montant ou descendant, appeie discontinuite dlnten- 
site, apparaissant sur la courbe dudit signal temporel a lisser. 

25 Dans le premier algorithms le gain de Kaiman est choisi completemerrt constant : il en resulte un trainage expo- 
nentiel affectant ledrt front de discontinuity dlntensite resultant d'un mouvement Ainsi, dans Ilmage d'origine bruitee, 
un petit objet, comme par exemple un catheter, qui aurait bouge rapidement donnant lieu a un creneau du signal 
d'intensite, peut avoir disparu de Pimage fBtree parce que les f lanes du creneau sont deformes par le filtrage. Cet algo- 
rithme efface les petits objets en mouvement. 

30 Dans le second algorithme, le gain de Kalman est une fonction de la difference entre I'intensite bruitee observee a 
un instant donne et I'intensite f iltree a I'instant precedent pour un pixel de localisation donnee. II en resulte que le signal 
temporel est lisse avant la discontinuite, mais que, par contre, il n'est plus fiitre apres la discontinuite ; il existe done un 
legs de bruit apres le front de discontinuite d'intensite. 

Le filtrage temporel connu a done le desavarrtage de ne pas pouvoir etre applique efficacement a une sequence 

35 dlmages fortement bruitees representant des petits objets animes. 

Ainsi, le filtrage temporel connu ne resout pas certains problemes importants qui se posent lorsque le filtrage tem- 
porel est applique a une sequence dlmages obtenues en mode f luoroscopique de rayons X, utiiisee, par exemple, pour 
suivre en temps reel une operation medicale durant laquelle un outii de diametre extremement petit, tei qu'un catheter, 
est introduit ou deplace dans la zone observee: 

40 

SUMMARY OF THE INVENTION 

Un probleme est que. du a la tres faibie intensite du faisceau fluoroscopique de rayons X, les images de la 
sequence sont extremement bruitees, et contiennent frequemment des pics de bruit. 
45 Un autre probleme est que, du au fait que chaque image de la sequence est separee de la suivante par un petit 
laps de temps, un evenement important comme le deplacement d'un petit outii tel qu'un catheter, peut survenir d'une 
image a I'autre. Le filtrage de Ilmage ou apparait le mouvement de ce petit objet ne doit pas deformer ou effacer cet 
objet. 

La presente invention vise a fournir un precede de filtrage temporel pour la reduction du bruit, dans les images sue- 
so cessives d'une sequence d'images numerisees : 

qui travailie strictement en temps reel, e'est-a-dire qui fournisse Timage courante filtree des I'apparition de ladite 
image courante observee, immediatement, sans retard dlmage, quelle que sort la cadence des images dans la 
sequence, 

qui reduise les legs de bruit suivant un front de discontinuite du signal temporel dlntensite, sans amortir le front 
55 de discontinuite ; 

qui sort capable de reduire les pics de bruit ; 
qui n'efface pas et ne deforme pas les petits objets en mouvement 
Ces buts sont atteirrts au moyen d'un procede de filtrage temporel du bruit dans une image appelee image du pre- 
sent d'une sequence d'images en forme de matrices bidimensionnelles de pixels, lesquels ont des vaieurs dlntensite, 
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appeiees echarrtillons bruitees numerisees, 

ce procede comprenant revaluation cfun echantillon filtre appele echantillon causa! filtre du present, pour 
reconstruire un echantillon bruite correspondant a un pixel de localisation donne, dans Timage du present par une 
combinaison lineaire de I'echantillon brurte du present et des echantiilons anterieurs dans la sequence, appeles echan- 
s tillons causals bruites affectes de poids qui sont des coefficients de probabilite de continuite dlntensite. de 0 (zero) a 
1; entre ledit echantiilon causal brurte auquei le poids est affecte et I'echantillon bruite du present le poids affecte a 
I'echantillon bruite du present ayant la valeur 1 . 

Un dispositif pour mettre en oeuvre ce precede de f iltrage comprend : 

10" un systeme de traitement dlmages pour fournir une intensite numerisee bruitee appelee echantillon du present 
brurte, d'un pixel ayant une localisation donnee dans une image ayant la forme d'une matrice de pixels arrivant a 
cet instant du present, 

et un sous-ensemble, appele sous-ensemble causal, recevant a son entree I'echantillon brurte du present, incluant 
une memoire pour fournir une intensite numerisee appelee echantillon causal bruite ou filtre du pixel anterieur au 
15 pixel du present, de meme localisation dans la matrice de Timage arrterieure. une memoire pour fournir un facteur 
de gain causal, et des moyens de calcul pour evaluer ('equation dlntegration causale, et pour fournir a sa sortie la 
valeur de I'echantillon causal filtre du present, qui est I'echantillon filtre pour reconstruire Techantiilon du present 
bruite. 

20 Ce dispositif off re I'avantage d'etre simple a mettre en oeuvre et de fournir strictement et exactement en temps reel 
un filtrage temporel du bruit efficace qui ne detruit pas les details de Timage. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 L'invention est decrite ci-apres en detail, en reference avec les figures schematiques annexees dont : 

la FIG.1 qui represente un dispositif de prises de vues a rayons X ; 
la F1G.2 qui represente une sequence dlmages numeriques ; 

la FIG.3A qui represente un signal temporel brurte correspondant a un pixel A localise aux coordonnees x t y, avec 
30 une discontinuity due a un mouvement ; 

la F1G.3B qui represente un autre signal temporel bruite, avec un pic de bruit ; 

!a FIG.4A qui represente un signal temporel filtre correspondant au signal temporel bruite de la FIG.3A ; 
la FIG.4B qui represente un signal temporel filtre correspondant au signal brurte de la FIG.3B ; 
la F1G.4C qui represente un signal temporel filtre correspondant au signal temporel bruite de la FIG.3A, selon une 
35 variante de ('invention ; 

les F1G.5A a 5C qui illustrent une determination de I'ecart type du bruit a B par rapport a la moyenne du bruit m B ; 
les FIG.6A a 6D qui montrent divers exemples possibles de la fonction F C (Z C ) ; 

la FIG.7 qui montre sous forme de blocs fonctionnels. un dispositif pour mettre en oeuvre le procede de filtrage 
tempore!. 

40 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
1/ Dispositif a rayons X 

45 En reference avec la FIG.1, un systeme de radiographie numerisee comprend une source de rayons X ; une table 
mobile 2 pour recevoir un patient ; un dispositif intensif icateur dlmage 3 couple a un tube video 4, qui fournit des don- 
nees a un systeme 5 de traitement numerique d'image comprenant un microprocesseur. Ce dernier comprend plu- 
sieurs sorties dont une sortie 6 couplee a un moniteur 7 pour la visualisation de la sequence dlmages radiographiques 
ou des images d'intensite. 

so L'image radiographique numerisee peut contenir 512 x 512 ou 1024 x 1024 pixels codes sur Sbits ou 10 bits. Cha- 
que pixel peut §tre ainsi affecte d'un parmi 256 ou 1024 niveaux d'intensite. Par exempie, les regions sombres ont un 
niveau d'intensite faibie, et les regions claires de Timage ont un niveau dlntensite eleve. 

Llmage numerisee peut etre obtenue en mode fluoroscopique. L'invention peut trouver son application particuiie- 
rement dans le filtrage d'images angiographiques. 

55 Dans tous les cas, la presente invention ne tient pas compte de la methode par laquelle la sequence dlmages 
numeriques a ete obtenue, ni de la nature des objets qu'elle represente, mais concerne seulement le filtrage de ces 
images dans cette sequence pour eliminer le bruit 
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It/ Precede de fOtrage tempore! du bruit 

La preserrte invention propose un procede de filtrage temporel du bruit successivemerrt de chaque image d'une 
sequence dlmages numeriques bnitees. Ce procede opfcre le filtrage de ta derni&e image bruitSe observee. Ce pro- 
5 ced6 est done exactement mene k bien en temps reel. , 

Ce procedS est capable d'effectuer ta detection de mouvement d'objets enregistres sur tes images de la sequence. 
II est fonde sur des etapes de filtrage rgcursif adaptatif. 

En reference avec la figure 2, le procede de ffltrage temporel du bruit comporte d'abord I'acquisition et la numeri- 
sation d'une sequence dlmages bruit ees, reaJisee a une cadence donnee. 
10 ' Cette sequence est composee, de llnstant le plus 6!oigne dans le passS. jusqu'& llnstant le plus recent : 

d'un nombre dlmages dlntensite du passe, notees Jj C appe!6es "images causales", formees depuis le premier 
instant j = t-k 0 , ou ko est un nombre entier egal au nombre dlmages de la sequence moins un, jusqu'a llnstant j = t-1 , 
ou t est llnstant ou arrive 11 mage en cours de filtrage ; 

d'une image J t p en cours de filtrage, appetee image du present, arrivant k llnstant j = t . 
is Chaque image numerisee bruitee Jj est une matrice bidimensionnelle de pixels not6s Aj (x, y), caracteris6 chacun 
par ses coordonnees x, y dans la matrice et par un niveau dlntensrte lj(x, y) sur une 6chelle de niveaux d'intensite, ou 
j est llndice correspondant k celui de llmage. 

Un signal, appel6 signal temporel bruite, est forme des differents niveaux dlntensite lj (x, y) des pixels Aj (x, y) en 
fonction du temps t (lire tau), comrne repr6sente sur les FIG.3A et 3B, ou I'intensite I est portee en ordonnee et x en 
20 abscisse. Le procede de filtrage temporel selon [Invention a pour objet le filtrage ou lissage du bruit qui affecte ce signal 
temporel note I(t), de mantere k obtenir un signal temporel filtre. note P(x), comme represents sur les FIG.4A et 4C. On 
appelle ci-apr6s "echantillons" ces intensity qui sont les points du signal temporel I(t). Les intensity du passe ou cau- 
sales sont notees lj C et tintensite du present est notee l p . 

Le procede selon llnvention s'applique au cas ou les grands mouvements se produisant dans la sequence d' image 
25 ont d6j& &6 compenses par des methodes de compensation de mouvement connus de I'etat de la technique, Par 
exemple des methodes de correction de translations. On evitera cependant les methodes d'estimation de mouvement 
couteuses en temps de calcul. Ce genre de methode doit etre evitee puisque rinvention vise la fourniture d'un procede 
applicable strictement en temps reel. 

I! faut noter que, m£me lorsqu'un procede de compensation de mouvement a ete ainsi applique sur une sequence 
30 d'images, on ne peut jamais affirmer que les mouvements aient ete reellement tous compenses, parce que les techni- 
ques de compensation de mouvement ne sont pas pariaites. 

Pour prendre en compte cet etat de fait, le present procede effectue un filtrage du bruit du signal temporel I(t) en 
traitant les petits mouvements locaux residuels ou non compensables des objets dans la sequence dlmage. 

Ce procede effectue un filtrage de I'echantillon bruite l t p arrive k Hnstant du present t, et construit un echantillon 
35 filtre P t c en utilisant I'echantillon bruite l t p du present t. et les echantillons anterieurs appeles echantillons causals brui- 
t6s l t -w) .... lt% . It% , Im . observes aux instants du passe de t-ko k t-1. 

Le filtrage temporel est effectue favorablement pour chaque pixel Aj (x, y) de la matrice bidimensionnelle, aux dif- 
ferentes coordonnees x, y, individuellement. 

La figure 3A montre un signal temporel bruite ayant, k titre d'exempie, une discontinuity dlntensite D entre 
40 I'echantillon 1 1 ?3 et I'echantillon 1 1 % due k un petit mouvement local. Ce petit mouvement se traduit par un front de mon- 
tee du signal temporel l(x) entre llnstant t-2 et llnstant t-3, un "haut plateau" de signal entre les instants t-3 et t, et une 
"plaine basse" aux instants precedant llnstant t-4. 

Le filtrage temporel selon llnvention est capable de s'adapter & une discontinue dlntensite telle que D. Ainsi, les 
F1G.4A et 4C morrtrent le signal temporel P(t) lisse par le procede selon llnvention dans diverses formes de realisation. 
45 Ce signal temporel montre en outre des petites dents de scie entre les instants t-8 et t-4, car, du au bruit, le niveau 
dlntensite d'un pixel de m§mes coordonnees (x,y) dans les images de la sequence varie constamment, provoquant un 
aspect de chute de neige. Le filtrage permet de supprimer cet aspect en lissant a une valeur moyenne constante le 
signal temporel l(x) relatif a un pixel donne. 

La FIG.3B montre par ailleurs un signal temporel bruite I(t) ayant a titre d'exemple un pic de bruit D* apparaissant 
so entre les instants t-4 et t-2, e'est-^-dire a I'instant t-3. Un pic de bruit tei que D* se distingue d'une discontinuite dlnten- 
site telle que D de la FIG.3A, par le fait qull ne correspond pas a un phenomfcne de mouvement spatial r6el. Un tel pic 
de bruit peut etre f iltr6 par le proc6d6 propose selon rinvention. 

L'Schantillon filtre P t c correspondant k I'echantillon bruite du present l, p s'exprime par la relation (1) : 

55 
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H-1 

Ib] C *'i C + '« P 

P . C = fe±!! Tn (1) 

I b i C + 1 

L'echantillon filtre P t c , ou signal temporel filtre a llnstant t, est une combinaison lineaire des echantillons bruites obser- 
ves aux instants causals de j = t-k 0 jusqu'a j = M , et de l'echantillon bruite observe a rinstant t. 

Dans la formuie (1) du signal filtre P t c . on affecte aux echantillons causals l£ , l t ? 2 etc.., un poids bf. 
Le denominateur de la formuie (1) est un facteur de normalisation dont la presence est fondee sur ie fait que la 
somme des poids que Ton applique aux differents echantillons doit §tre egaie a 1 pour que la moyenne du signal filtre 
P t c soit egale a la moyenne du signal bruite l t p . 

Le precede de fittrage temporel est fonde sur une combinaison lineaire des echantillons causals, avec rechantiilon 
is du present, dans laquelle les poids bj C , affectes aux echantillons causals ont une forme specifique, le poids affecte a 
l'echantillon du present etant 1 . 

Le poids affecte a un echantillon causal donne est une probability de continuity d'intenstte entre ledit echantillon 
causal donne et l'echantillon du present. 

D'une maniere generale, la formulation du poids causal b t c k a affecter a un echantillon causal l t c k estdonnee par 
20 la relation (2) : 


10 


b t°k = ft a< i ' aveck >1 (2) 

t-K+1 

25 

ouaj est la probabiiite de continuity d'intensite entre les echantillons causals successifs depuis l t Sc jusqu'a l t p . La 
relation (2) exprime qu'un poids b x % relatrf a un echantillon causal I & est egal au produit de toutes les probabiiites de 
continuity dlntensity a f depuis rinstant causal j=t-k+1 jusqu'a I'instant du present t. 

30 Ainsi, la determination du poids affecte a un echantillon causal donne de la formulation (1) est fondee sur I'hypo- 
these que ledit echantillon donne est retenu et pris en compte dans la combinaison lineaire seulement dans la mesure 
ou ii n'est pas trop different de l'echantillon du present ceia signifiarrt quil apparent au meme objet 

Selon cette hypothese, par exemple, la formulation du poids causal b t ?i , relatif a rechantiilon causal l t ?i , se fait 
en ecrivant que b t ?i est une fonction de la valeur absolue de la difference entre rechantiilon l t p a llnstant present et 

35 l'echantillon I m a llnstant precedent. 

De preference, le poids causal bt?i est une fonction de la difference : 


A c = |l t p - P^l (6 C ) 

40 

ou Pm est l'echantillon deja filtre a I'instant precedent t-1. et done repute moins bruite que l t o du faitdu filtrage. Si la 
difference entre les echantillons l t p et P t ?i estfaible, on affecte au poids correspondant b t?i une valeur "elevee" proche 
de 1 . Si cette difference est grande. alors on affecte au poids b t ?i une valeur proche de 0 (zero). Dans ce cas, I'echan- 

45 tillon I t ?i n'est pratiquement pas pris en compte. 

Ensuite, dans cet exemple, la formulation du second poids bt?2 causal relatif a l'echantillon causal l t ?2 se fait en 
ecrivant que le poids bt?a est une fonction non seulement de la difference entre rechantiilon a i'instant t et l'echantillon 
a llnstant t-1 . mais aussi est une fonction de la difference entre l'echantillon a I'instant t-1 . et l'echantillon a I'instant t-2. 
Ainsi, le poids b t ?2 a une valeur elevee proche de 1 , si on ne prend en compte que des echantillons n'ayant pas subi 

so de modifications autres que celles dues au bruit par rapport a l'echantillon du present l p . e'est-a-dire si I'on pose 
comme condition que les differences entre les echantillons pris en compte sont faibles. Ceci conduit a la formulation 
des poids causals bj C comme des produits de fonctions des cfifferences d'intensite apparaissant sur le signal temporel, 
e'est-a-dire a la formulation de ces poids comme des produits de probabiiite de continuity d'intensite entre les echan- 
tillons precedant l'echantillon du present a filtrer. 

55 Ainsi, la formulation du poids causal relatif a l'echantillon l t ?i est : 
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ht-x = «t (2a) 

oua? est une probability de continuity d'intensity entre rychantillon du present l t p et rychantillon ftltre P t ?i . La conti- 
nuity cfirrtensity entre rychantillon k llnstant t et rychantillon a llnstant pr6c6dent t-1 se traduit par une probability de 
continuity forte. Et une discontinuity dlntensite entre les echantillons se traduit par a? proche de 0. 
La formulation du poids causal relatif ^rechantillon ! t ?i est: 

b t c . 2 = a£ x (2b) 


is ou a t ?i est la probability de continuity d'intensity entre I'echantillon a llnstant t-1 et rychantillon a llnstant t-2. 

La continuity d'intensity des echantillons entre llnstant t et llnstant t-2 implique la continuity dlntensity entre les 
instants t et t-1 , et la continuity d'intensity entre les instants t-1 et t-2, ce qui se traduit par fe produit de probabilitys (2b) 
dans lequel la probability de continuity d'intensite entre les instants t et t-1 est forte, et la probability de continuity 
d'intensite entre les instants t-1 et t-2 est forte. 

20 La formulation du poids causal relatif a rychantillon l t % est ensuite : 

b t C -3 = ttt- 2 * at-! x a t (2c) etc. 

25 Ainsi dans I'exemple de la F1G.3A, du fait de la discontinuity dlrrtensite D entre les echantillons I & et \x% . la pro- 
bability de continuity a t % est tres proche de 0 (zero). Le poids b A obtenu par le produit des probabilites incluant cette 
probability <z t S tend ygalement vers zero. Ensuite le poids b x ? 5 tend aussi vers zero du fart que sa formulation est un 
produit contenant aussi cette probability proche de zero. Des qu'il y a une probability proche de zero dans le produit 
des probabilitys fournissant un poids bj C , la presence de cette probability annule ce produit et tous les produrts corres- 

30 pondarrt aux instants qui precedent I'instant ou cette probability proche de zero apparait II en resulte que les poids bj C 
des echantillons causals correspondant auxdits instants precydents sont aussi proches de zero. Dans I'exemple parti- 
cuiier de la FIG. 3A, les poids causals sont : b t ?i - 1 ; b t ?2 - 1 ; b t S - 1 ; b& - 0 ; b t % - 0 

If en resulte que, par cette formulation spycifique (2) des coefficients causals bj C , ryquation (1), dont le calcui etait 
tres complexe, est transformee en un seul groupe de 3 relations simples realisant le firtrage temporel recherche. 

35 L'yquation lineaire (1), dans laquelle les poids causals ont une forme specifique (2) conduit a une formulation de 
Tequation recursive (3 C ) suivante : 

P C C = P'-x + K t C X (I t P - P^) (3 C ) 

40 

Cette equation (3 C ) appelee yquation dlntegration, prend en compte rechantillon causal filtry P t ?i a I'instant t-1, 
modrfie par ce qu'on appelle ci-apres un terme d'innovation causal constituy par la difference entre rychantillon filtre 
P t ?i a i'instant t-1 et rychantillon bruite l t p observe a I'instant du present t Ce terme d'innovation est multiplie par un 
terme appele gain causal K t c , qui varie entre 0 et 1. Cette prendre equation (3 C ) fournit une valeur filtree P t c qui est 
45 I'echantillon filtry recherche. 

Le gain causal K t c n'est pas un gain de Kalman, mais un facteur de gain completement diffyrent donne par Tequa- 
tion recursive suivante (4 C ) : 

K t C -i + *t 

Dans le passage entre les equations (1 ) et (2) et les equations (3 C ) et (4 C ), le gain causal est Knverse de la somme 
55 des poids causals 
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K t C = 
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K t c = 1/ J bj C 

j=t-Ko 

5 > 9 * I HI 

et de plus, comme on vu plus haut : 

*t - b t c . x (2a) 

10 

qui exprime que la probability de continuity a? est le poids de I'echantiHon causal I t ?i precedant Pecharrtillon l t p du 
present probability qui est appelee ci-apres coefficient causal de continuity. 

Le passage entre les equations (1) + (2) et I'Squation (3 C ) est alors fonde sur le systeme d'equations causales : 

15 P t c = Pt 0 -! + K t c x (I t - P^) equation d 1 integration (3°) 

K c c = — facteur de gain causal (4^) 

Kt C -i + «t 


et = F c ( ) coefficient causal de continuite (5 ) 

s t c 

25 

ou S t c est un facteur de normalisation, et ou F° est une fonction definie plus loin. 

^iteration de I'equation (3 C ). en commencant a t-ko dans le temps, permet de retrouver exactement I'equation (1 ), 
telle que formulee plus haut 

30 La formulation (3 C ) de I'equation d'integration correspond done directement a cette combinaison lineaire (1) des 
echantillons associes a des poids specifiques qui sont des fonctions de continuite d'intensite entre les echantillons sue- 
cessifs et qui sont determines par le produit des probabilites de continuite (2) entre les echantillons successrfs, de 
I'echantilton dont fe poids est calcui e jusqu'a I'^chantillon du present. 

35 Formulation du coefficient de discontinuity causal a? 

Le coefficient causal de continuite a? est defini comme une fonction de la difference 

A c = |I C P - P'-il (6 C ). 

I! convient de determiner si la difference A c est relative a une discontinuite telle que D sur la FIG.3A, ou si elle est 

seulement relative a du bruit 

Pour determiner la participation du bruit, la difference A c est normalisee par un facteur S t c prenant en compte la 
45 variance du bruit relative a chaque echantillon de cette discontinuite A c . L'ecart type du bruit note a B , mesure en 

niveaux d'intensite, pour chaque echantilion, peut etre estime par toute methode connue de Tetat de la technique, ou 

bien peut §tre estime a priori. 

On rappelle que la variance a| du bruit est la moyenne des ecarts au carre du bruit dans une region consideree 

de Timage. les ecart dlntensite dus au bmit etant mesures par rapport a la moyenne d'intensite locale. 
so La determination de l'ecart type du bruit peut se trouver compiiquee du fait que dans le systeme de la FIG.1 on 

introduit une transformation de llmage selon une loi de compression des niveaux d'intensite faibles et eleves avant la 

numerisation du signal pour tenir compte du fait que le bruit est moins visible a I'oeil dans les regions tres sombres ou 

tres brillantes, que dans les regions d'intensite moyenne. L'ecart type du bruit a B est alors une fonction de llntensite a 

chaque instant, en chaque pixel. 
55 La determination de l'ecart type du bruit doit done etre faite localement en chacun des pixels de llmage, en fonction 

de la moyenne d'intensite de la region correspondante : ainsi le bruit est caracterise spatialement du fait que Ton calcule 

localement une variance du bruit et une moyenne du signal observe, et on arrive done a separer ce qui est du a llmage 

meme, ou signal, du bruit 

Selon une methode proposee a titre d'exemple, en reference avec la FIG.5A, on determine dans une image Jj des 
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petites fenStre Wj, typiquement des zones de 5 x 5 ou.3 x 3 pixels autour de chaque pixel Aj(x.y). Et dans chacune de 
ces petites fenStres Wj(x,y), on determine la variance du bruit et la moyenne d'intensite observe. 

A cet effet n etant le nombre de pixels d'une petite fenetre Wj(x,y), pour chaque position de cette petite fen§tre cen- 
tree sur un pixel Aj(x,y), on determine : 

les intertsites lj(KO de chacun des pixels de la petite fenetre Wj(x,y), les indices k et / etant les coordonnees des 
pixels considered situes dans cette fenStre par rapport aux coordonnees rectangulaires X et Y de Hmage Jj ; 

la moyenne d'intensite my^x.y) dans cette fenelre, en cette position, donnee par la relation : 


la variance du bruit a BW {x,y) au carre, dans cette fenetre, donnee par la relation : 
15 <4,(x.y) 2 = ± E [I.(k.2) - m w (x.y)] 2 (9°) 


ou interviennent sous le signe somme L la difference au carre entre les intensites de chaque pixel de la fenetre et la 
20 moyenne rr^(x t y) ; on prend en outre la moyenne de cette somme par le nombre n de pixels dans la fenStre. 

En reference avec la FIG.5B, on realise ensuite un histogramme a deux dimensions dans iequel on porte sur I'axe 

horizontal des valeurs moyennes mw(xy) de Tintensite avec une precision donnee : par exemple de 10 niveaux de gris 

en 10 niveaux de gris, et sur I'axe vertical la valeur moyenne de la variance. Dans chaque intervalle de la moyenne 

d'intensite mw(x,y) porte sur I'axe horizontal, on trouve le decompte de la frequence N de la variance ou nombre de fois 
25 qu'une variance donnee se produit pour cet intervalle de la moyenne dlntensite. Tout a fait a titre d'exemple, pour la 

ciarte de I'expose, sur la FIG.7B, on a represents qu'entre les valeurs 10 et 20 de la moyenne m w (x,y): 

la frequence N 01 de la variance ayant la valeur 1 est 0 

la frequence N 1 1 de la variance ayant la valeur 4 est 10 

la frequence N 2 i de la variance ayant la valeur 9 est 20 
30 la frequence N 31 de la variance ayant la valeur 16 est 35 

la frequence N 41 de la variance ayant la valeur 25 est 

etc.. 

La methode proposee se poursuit en reference avec la FIG.5C par le trace de la courbe de la frequence N d'occur- 
rence des variances portee sur I'axe vertical en fonction de la variance a|w (x,y) , pour tout intervalle considere de 
35 Pintensite moyenne m^y). Cette courbe des frequences N en fonction de la variance passe par un maximum qui 
donne une valeur possible de la variance, et qui doit etre determine pour chacun des interval les de la moyenne d'inten- 
site. 

On obtient ainsi une variance du bruit en fonction de la valeur moyenne d'intensite attribuee a chaque pixel Aj(x,y). 
Cette technique permetd'extraire rinformation concernant le bruit a partir de I'image m§me. 
40 L'ecart type du bruit a B est alors obtenu par ie calcul de la racine carree de la variance du bruit : 


j=t-N+l 


La variance de la contribution du bruit entrant dans la difference A c peut §tre formulee par : 


2-K t C , 


en regime stationnaire ; cette valeur est utilisee comme une approximation pour la normalisation de A c . Le facteur de 
normalisation de la difference A c est alors donnee par la relation : 


55 
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S t c = cr B —L^ (11*) 

Dans r equation (5°). le coefficient de continuite a? est done une fonction d'une difference & c donnee par la rela- 
tion (6°), normalisee par une approximation St C de l'ecart type du bruit <j Bi pour prendre en compte le fait que la diffe- 
rence A c est affectee par le bruit A cet effet, le caicui du coefficient a? comprend, en chaque point de la courbe I(t) 

io du signal tempore!, la determination de l'ecart du bruit a B par rapport k la moyenne. 

La force du procedede fUtrage propose selon I'invention, fonde sur le caicui de la fonction d'integration causale (3 C ) 
associee aux relations (4 C ) et (5°). reside tout particulierement dans le fait que Texpression a t est normalisee comme 
propose. Puisque Ton a pris en compte, pour cette normalisation non seuiement la variance du bruit de I'echantillon du 
present l t p mais aussi la variance du bruit de I'echantillon f iltre P t ?i , I'equation (5°) donnant le facteur de gain, per met 

is une estimation du bruit residuel pour chaque echantillon f iltre. 

La fonction d'integration recursive (3 C ) est maintenant aisement calculable car elle depend uniquement du facteur 
de gain K t c donne par I'equation (4 C ). Ce facteur de gain K t c depend lui-meme uniquement du facteur de gain K t ?i deja 
calcule dans le passe et du coefficient de continuite a? . Ce coefficient de continuity a? depend aussi du facteur de 
gain deja calcule Km , d'ou il resutte que ce coefficient de continuite a? doit etre d'abord calcule, que I'equation de 

20 gain K t c est ensuite aisement calculable avec le coefficient de continuite a? , et que I'echantillon f iltre P t c est ensuite 
aisement calculable par ia relation (3 C ). 


Definition de la fonction F° 


. 25 Pour le caicui de i'equation (5 C ) donnant a? , la fonction F° est maintenant definie ci-apres. 
L'argument de la fonction F° est note z c : 


z c = 


I P -P c 


t-i 


(12 c ) 


30 


35 

La fonction F C (2 C ) est tracee a titre d'exemple sur les FIG.6A k 6D en portant ses valeurs en ordonnee, et en por- 
tant z c en abscisse. 

La fonction F(z°) a d'abord une valeur constante Fornax lorsque 0 * z c ^ 1 ; puis la fonction F^z 0 ) est decroissante 
jusqu'& une valeur r^min pour 1 < z c . La valeur Fmax est telle que : r^max ^ 1. La valeur r^min est telle que 0 ^ 
40 F c min. Des exemples de teiles fonctions F(z°) sont donnes ci-apres. 

Exemple I : illustre par la RG.6A 

La fonction Ffz 0 ) est constante et egale k une valeur F c max = 1 quand 0 * z° * 1 , puis la fonction F c (z c ) est 
45 choisie decroissante jusqu'a une valeur F c min = 0 quand 1 s z c . La fonction F^z 0 ) est avantageusement, dans cette 
region, de la forme d'une Gaussienne affectee d'un facteur de normalisation, telle que : quand 1 £ z c . alors 



50 

Exemple II : illustre par la FIG.6B 

La fonction F c (z c ) est constante et egale k une valeur F c max inf erieure a 1 , par exemple F c max = 0,85 quand 0 £ 
z c < 1 ; puis la fonction Ffz 0 ) est decroissante quand 1 s z c . La fonction F c (z c ) est limitee en decroissance par exemple 
55 k = 0, 1 0. La fonction F^z*) peut £tre, dans ce cas, une Gaussienne decalee. 

Exemple HI : illustre par la F1G.6C 

La fonction F c (z c ) est constante et egale k une valeur Fornax inferieure a 1, mais proche de 1. telle que F c max = 


9 


EP0751 482 A1 


0,85 lorsque 0§z c 5l. Puis la fonction F^z) est lineairement decroissante jusqu'a une valeur Fmin proche de 0, par 
exemple ^min = 0,10, qu'eile atteint par exemple pour z c = 2. Ensuite la fonction F^z 0 ) est lineairement cortstante et 
egalea 

c C 

5 i , • i F min = 0,10 pour 2 <z 

Exemole IV : illustre par la FIG.6D 

„ La fonction F c (z c ) est une fonction lineaire par morceaux qui est une approximation de la fonction de rexemple I, 
10 ou bien de la fonction de rexemple II. 

La fonction F^z 0 ) peut encore §tre choisie parmi cfautres formes appropriees dont la definition est a la portee de 
rhomme du metier. 

La fonction F^z 0 ) regie la puissance de filtrage. Lorsque F^nnax = 1 . la puissance de f iltrage est maximaie Lorsque 
Fornax est inferieure a 1, la puissance de filtrage est legerement limited Lorsque F°min est superieure a zero, la puis- 
15 sance de filtrage n'est jamais nulle. 

Exemple de realisation de llnvention avec une fonction F° de rexemple I 

Les FIG.4A et 4C montrent chacune un signal tempore! filtre P(x) obtenu a partir d'un signal tempore! bruite I(t) 
20 ayant au moins une discontinuity d'intensite D due a un mouvement, tel qulllustre par la FIG.3A. 

Ce filtrage est effectue au moyen du systeme de filtrage (3 C ), (4 C ). (5°) dans les conditions ou F°(2?) est une fonc- 
tion analogue a celle de I'exempie I illustree par la FIG.6A. 

En reference avec la FIG.4A, le filtrage causal produit un lissage du signal tempore!, a gauche de la discontinuity 
D, dans ie sens causal de t-4 vers t-ko ; ensuite, du fait de la discontinuity D, le filtrage causa! n'est pas capable de rea- 
25 liser le filtrage de rechantillon bruite \ x % en t-3, car il a "oublie" le passe constitue par les echantiilons l t ^ , It-s etc ; 
de sorte que le filtrage causal est tres efficace en de la courbe de signal temporel et a gauche de cet echantillon, 
puis le signal est directement le signal bruite 1 1 % , en t-3. car en ce point, le filtrage causal ne fournit que I'intensite brui- 
• tee observee, non modrfiee par les donnees du passe. Ensuite ie filtrage causal produit a nouveau un signal tres bien 
lisse a droite de I'echantillon l t % , a partir du point suivant, en t-2, du signal temporel. ou I'echantillon en t-3 peut etre 
30 pris en compte du fait qui I n'y a pas de discontinuity dlntensrte errtre les instants t-3 et t-2. ^ 

Done, il reste apres filtrage causal, tel que montre sur la FIG.4A, la premiere dent de bruit en l t -3 a I'instant t-3 du 
fait que les coefficients b£ affectys au signal en t-4 sont nuls puisque la difference d'intensite errtre t-3 et t-4 est 
grande. 

D'autre part, en refyrence avec la FIG.3B, le signal temporel bruite l(-c) peut montrer aussi un pic de bruit D' par 
35 exemple apparaissant a I'instant t-3. Le filtrage causal ne permet pas de distinguer la discontinuity D* due a un pic de 
bruit dans le signal temporel bruite de la FIG.3B, de la discontinuity D due a un mouvement dans le signal temporel 
bruite de la F1G.3A. Done le filtrage causal ne permet pas de bien filtrer un pic de bruit D' tel que montre sur la FIG.3B 
si la fonction choisie F C (Z C ) est celle de I'exempie I iilustrees par les FIG.7A. 

40 Exemple de realisation de rinvention avec une fonction F° de I'exempie H 

Selon I'invention, on peut modifier a volonte la puissance de filtrage en modifiant les valeurs de F c min, et F c max, 
de maniere a diminuer au mieux la dent dintensite bruitee residuelle, obtenue a droite de la discontinuity d'intensity D 
du signal temporel dans Texemple precedent tout en conservant un lissage aussi parfait que possible de ce signal tem- 
45 porel. 

Ainsi, selon rinvention, le passe peut n'etre oublie que "dans une certaine mesure", d*ou il resulte que le filtrage 
causal peut n'itre pas completement nul en l t ?3 - 

En reference avec la FIG.4C qui est le resultat du filtrage du signal de la FIG.3A en utilisant une fonction F c (z c ) 
ayant Fornax = 0,85 et F c min = 0,10, iilustrees par la F1G.6B, on montre que la premiere dent de bruit en l t S apres la 
so discontinuity D est amortie du fait que la fonction F 0 ^) n'est jamais nulle, puisque dans cet exemple au moins egale 
a 0,10. Done les donnees au temps t-4 sont toujours prises en compte et permettent de filtrer toujours dans une cer- 
taine mesure le signal au temps t-3, apres la discontinuity d'intensite D. 

Par ailleurs, comme on a vu precedemment, on peut aussi modifier a volonte la puissance de filtrage en modifiant 
les valeurs F° m in et F 0 ^ pour filtrer un pic de bruit D* tel que montre sur la FIG.3B. 
55 Ainsi, considerant la FIG.4B qui est le resultat du filtrage du signal de la FIG.3B en utilisant une fonction F C (Z C ) 
ayant F c max = 0,85 et F c min = 0.10, iilustrees par la F1G.6B, le pic de bruit D' est amorti du fait que la fonction F C (Z C ) 
n'est jamais nulle, car les donnees au temps t-4 sont toujours prises en compte et permettent de filtrer toujours dans 
. une certaine mesure le signal du pic de bruit D\ au temps t-3. 

On notera que la premiere image flltree est en fait la deuxieme image am vant dans la sequence puisque Ton a 
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besoin des donnees de ilmage precedente pour filtrer Ilmage courante. Les donnees de fa premiere image de la 
sequence ne sont pas filtrees pour leurs propres affichages et ne sort pas filtrees pour ('utilisation du fiftrage de Ilmage 
arrivant en deuxieme. Le filtrage reef commence a cette deuxieme image. 

5 III/ Dispositif pour mettre en .oeuvre le proced6 * 

En reference avec la FIG. 7, un dispositif simple pour mettre en oeuvre le precede de filtrage temporel precedem- 
ment decrit, est represents sous forme de blocs fonctionnels. 

Pour I'image du present a (Instant t Hntensite bruitee l t p d'un pixel de localisation x.y est amenee a 1'errtree 1 00 du 
to dispositif. 

Cette intensite bruitee l t p suit deux voies. 

Dans une voie, le signal bruite l t p est amene vers une fonction a B tabulae dans une Table LUT 1 1 (en anglais, LUT 
= LOOK UP TABLE) qui fournit i'ecart type du bruit c B du signal bruite l t p . L'ecart type du bruit a B est disponible a la 
sortie de la Table 1 1 . 

is Les gains K sont stockes dans la memoire MEM 17. Cette memoire MEM 17 a une sortie pour fournir les gains 
K t ?i de llnstant precedent, pour le pixel de localisation donnee x,y, et a une entree pour recevoir les gains K t c calcules 
pour I'instant du present t 

Le dispositif comprend une Table LUT 13 (LOOK UP TABLE) qui contient des vaieurs tabutees pour fournir imme- 
diatement la valeur 
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des qu'on lui introduit la valeur de a B venant de la Table LUT 1 1 , et la valeur du gain K t ?i venant de la memoire MEM 
17. La valeur constitue le denominateur de I'expression (5°) qui permet le caicul du coefficient de continuite a? . 
Cette valeur de denominateur S t c est disponible a la sortie de la Table LUT 13 et est introduite dans la Table LUT 14 
30 (LOOK UP TABLE). 

Le dispositif comprend une seconde memoire MEM 10 pour stacker ['intensity fittree en dernier pour le pixel traite. 

A llnstant du present t ou Ton cherche a calculer i'intensite filtree P t c , la memoire MEM 10 contient done I'intensite 
f iltree P ^ correspondant a Tinstant precedent et fournit done P t?i sur sa sortie. 

L'intensite f iftree a Tinstant precedent P t ?i est amenee sur I'additionneur 1 2 en meme temps que I'intensite bruite 
35 l t p observee a I'instant t. L'additionneur 12 fournit a sa sortie la difference entre les deux intensites Tune bruitee a I'ins- 
tant t notee ! t p Tautre filtree a I'instant t-1 , notee P t ?i . 

La Table LUT 14 prend la valeur absolue de la difference fournie par l'additionneur 12. et recoit en outre le deno- 
minateur S t c . La Table LUT 14 fournit alors, a sa sortie, I'argument z c de la fonction F°(Z) C : 


z c = Af . | J t p t°-i | (12 c } 


La Table LUT 15 (LOOK UP TABLE) contient la fonction F C (Z) C choisie tabulee. Cette Table LUT 15 recoit a son 
entree I'argument z° et fournit le coefficient de continuite a T 

s t 

La Table LUT 16 contient la formulation tabulee pour fournir ie gain a I'instant t a partir du gain de I'instant precedent. 


provenant de la memoire MEM 1 7, et a partir du coefficient de continuite a ? arrivant de la Table LUT 1 5. Ce gain 
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est fburni en sortie de la Table LUT 16, et est dirige d'une part vers la memoire MEM 17 pour remplacer ia valeur du 
gain K t ?i de Pinstant precedent, et d'autre part vers un muitiplicateur 19. 
■ La fonction d'integration (3°) peut maintenant §tre calcuiee. 
w L'additionneur 18 report d'abord le signal bruite observe a Hnstant t note l t p provenant de Tentree 100 sur une 
seconde voie, et le signal filtre k I'instant precedant t- 1 provenant de la sortie de la memoire MEM 1 0, en fournit la dif- 
ference l t p - Ph. Cette difference provenant de l'additionneur 18 est amenee sur le muitiplicateur 19 en nrteme temps 
que le gain K t c provenant de la Table LUT 16. Le muitiplicateur 19 effectue le produit 

is K^xO^-P^ ) 

et en fournit le resuftat a sa sortie. 

L'additionneur 20 regoit d'une part le produit venarrt du muitiplicateur 19, et cfautre part le signal deja filtre k I'ins- 
tant precedent, note P t ?i . provenant de la memoire MEM 10. L'additionneur fournit a la sortie generate 200 du filtre, la 
20 valeur de I'intensite filtree a I'instant t. Cette valeur filtree P, c vient aussi remplacer, dans la memoire MEM 10, llntensite 
firtree a I'instant M, qui est introduite par I'entree 120 de cette memoire 21 , selon I'equation d'integration : 

P t C = Pt C -i - K c c x (l t P - P^) (3^) 

25 

On notera que le disposrtif decrit plus haut, ne necessite la memorisation que de I'echantillon filtre P m et du gain K t ?i 
de I'instant precedent et ne necessite en outre que la valeur de I'intensite bruitee l t p a I'instant du present, pour mettre 
en oeuvre le precede de filtrage temporel. II n'est pas indispensable de memoriser llntensrte arrivant a I'entree 100. 
Apres filtrage des difterents pixels de I'image du present cette image est par exemple rendue visible sur le moni- 
30 teur 7 du disposrtif d'affichage de la FIG. 1 . Les differerrtes images de la sequence peuvent ainsi etre aff icrtees apr£s 
leur traitement respectrf. 

Revendications 


as 1. Procede de filtrage temporel du bruit dans une image (J^ appelee image du present, d'une sequence dlmages 
en forme de matrices bidimensionnelles de pixels, lesquels ont des valeurs d'intensite, appelees echantillons 
(lj(x,y)), bruitees numerisees, 

ce procede comprenant revaluation d'un echantillon filtre (P t c ) appele echantillon causal filtre du present, 
pour reconstruire un echantillon bruite (I,*) correspondant k un pixel de localisation (x,y) donne, dans I'image du 

40 present par une combinaison lineaire de I'echantillon bruite du present (l t p ) et des echantillons anterieurs dans la 
sequence, appeies echantillons causals bruites (lj C ) affectes de poids (bj C ) qui sont des coefficients de probabilite 
de continuite d'intensite, de 0 (zero) a 1, entre ledit echantillon causal bruite auquel le poids est affecte et I'echan- 
tillon bruite du present, le poids affecte k I'echantillon bruite du present ayant la valeur 1. 

45 2. Precede selon la revendication 1 , selon lequel le poids affecte k un echantillon bruite causal donne est le produit 
•des coefficients de probabilite de continuite entre les echantillons causals bruites successifs, de I'echantillon bruite 
causal donne jusqu'& I'echantillon bruite du present, et precede selon lequel la combinaison lineaire des Echan- 
tillons bruites causals et du present est normalisee par la somme des poids relatifs auxdits Echantillons. 

so 3. Procede selon la revendication 2, selon lequel. pour ('evaluation de i'echantillon causal filtre (P t c ) du present (t), la 
combinaison lineaire est realisee en effectuant revaluation d'une relation recursive, appetee relation d'integration 
causale, comprenant la determination d'une somme de TEchantillon causal filtre (Pm ) de (Instant precedent (t-1) 
et du produit d'un facteur, appele facteur de gain causal (K t c ) , par une difference entre I'echantillon bruite du pre- 
sent (ift et I'echantillon causal filtre (P t ?i ) de I'instant precedent 

55 

4. Procede selon la revendication 3 selon lequel. ledit facteur (K t c ) de gain causal de I'instant du present (t) est evalue 
de maniere recursive, par le rapport du facteur de gain causal (K t ?i ) de llnstant precedent (t-1) sur la somme ducfit 
facteur de gain causal (K m ) de I'instant precedent et d'un coefficient appele coefficient de continuite causal (a t ). 


12 


EP0 751 482 A1 


5. Precede selon la revendication 4, selon lequei ledit facteur de gain causal (K,) est Hnverse de la somme des poids 
affect es aux echantilions causais de la combinaison lineaire, et le coefficient de continuite causal (a t ) est le poids 
de recharrtillon causal (Im ) precedent rechantillon du present (O- 

5 6- Precede selon Tune des revendications 4 ou 5, selon lequei revaluation de la relation dlntegration causale com- 
prend : 

la determination de la valeur absolue d'une difference d'echantillons appelee difference causale (A c ) entre 
rechantillon bruite de Hnstant du present (t) et rechantillon causal fittre (P m ) de (Instant precedent (M ), 
w et la determination du coefficient causal de continuite (a ? ) comme une fonction decroissante (F°) de ladrte 

difference causale (A c ). 

7. Precede selon la revendication 6, seion lequei. pour revaluation du coefficient de continuite causal (a? ) ? la fonc- 
tion d6croissante est une fonction appelee causale (F^) dont rargument (Z 0 ) est donne par la difference causale 

is (dp), normalis6e par un facteur (S, c ) egal a la racine carree du cumul des variances du bruit de chacun des echan- 
tiilons de la difference, par rapport a la moyenne du bruit 

8. Precede selon la revendication 7, seion lequei la fonction decroissante causale (F°) a une valeur maximale cons- 
tante inferieure ou egale a 1 lorsque son argument (z c ) est compris entre 0 a 1 , et decroit vers une valeur minimale 

20 egale ou sup6rieure a zero, iorsque son argument est superieur a 1 . 

9. Precede selon la revendication 8. selon lequei la puissance de fittrage est regiee par les valeurs maximale et mini- 
male de la fonction decroissante causale (f°) choisies pour lisser les pics de bruit et pour lisser les dents de bruit 
qui suivent les discontinues dues a un mouvement spatial. 
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10. Dispositif pour mettre en oeuvre un precede de fiftrage selon Tune des revendications 1 a 9, comprenant : 

un systeme de traitement dlmages (5) pour fournir une intensite numerisee bruitee (1^ appelee echantiilon du 
present bruite, d'un pixel [A t (x,y)] ayant une localisation donnee (x.y) dans une image (J^ ayant la forme d'une 

30 matrice de pixels arrtvarrt a cet instant du present (t), 

et un sous-ensemble, appeie sous-ensemble causal, recevant a son entree (100) recharrtillon bruite du pre- 
sent (l/ 3 ), incluant une memoire (10) pour fournir une intensite numerisee appelee echantiilon causal (l t ?i ; 
P t ?i ) bruite ou filtre du pixel anterieur au pixel du present, de m£me localisation (x.y) dans la matrice de 
llmage arrterieure, une memoire (17) pour fournir un facteur de gain causal (Km ). et des moyens de calcul 

35 pour evaluer requation dlntegration causale. et pour fournir a sa sortie (200) la valeur de recharrtillon causal 

filtre (P t c ) du present, qui est l'6chantiilon filtre pour reconstruire recharrtillon du present bruite (I,**). 

11. Dispositif selon la revendication 10, dans lequei les moyens de calcul, pour revaluation de requation dlntegration 
causale et la fburniture de recharrtillon causal filtre (P t °) comprennent : 

40 

un additionneur (12) pour la determination d'une valeur absolue d'une difference d'echantillons appelee diffe- 
rence causale (A c ) entre recharrtillon bruite de rinstant du present (t) et recharrtillon causal filtre (Pm ) de 
I'instant precedent (t-1) memorise, et des tables (LUT 13. 14, 15) pour la determination d'un facteur de conti- 
nuite causal (a ? ) comme une fonction decroissante (F°) de la difference causale (A G ) normalisee par un fac- 
45 teur (S t c ) de la variance du bruit (o^), 

une table (LUT 16) pour la determination d'un facteur (K t c ) appeie facteur de gain causal du present (t) egal 
au quotient du facteur de gain causal (K t ?i ) de I'instant precedent (t-1) memorise, par la somme dudit facteur 
de gain causal (Km ) et du coefficient de discorrtinuite causal (a? ), 

des additionneurs (18, 20) et multiplicateur (19) pour la determination de la somme de rechantiilon causal filtr6 
so (P t ?! ) de I'instant precedent (M) et du produit dudit facteur de gain causal (K t c ) par une difference entre 

recharrtillon du present bruite (l t ^) et recharrtillon causal filtre (Pm ) de I'instant precedent, donnant ainsi la 
valeur de recharrtillon causal filtre (P, c ). 

12. * Dispositif pour I'aide a ('observation d'une image medicale, comprenant : 

55 

un systeme pour fournir des donnees d'une image sous forme de matrice bidimensionnelle de pixels ayant des 

valeurs dlntensite numerisees, 

un systeme d'affichage pour i'affichage de llmage, 

un systeme de traitement dlmages ayant acces aux donnees d'image, et au systeme daffichage, 


13 


w EP 0 751 482 A1 ^ 

ce syst&me de traitement dlmages incluant un dispositif seton I'une des revendications 10 ou 11 . 
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